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Abstract: A single impurity in the ID Luttinger model creates a local modification of 
the charge density analogous to the Friedel oscillations. In this paper, we present an exact 
solution of the case g = \ (the equivalent of the Toulouse point) at any temperature T and 
impurity coupling, expressing the charge density in terms of a hypergeometric function. 
We find in particular that at T — 0, the oscillatory part of the density goes as lnx at small 
distance and x -1 / 2 at large distance. 



06/96 



1. Introduction 



The Luttinger model, describing the low-energy excitations of an interacting one- 
dimensional fermion gas, is one of the simplest non-Fermi-liquid metals. Experimental 
observations of this non-Fermi state in ID quantum wires are difficult, since disorder tends 
to localize the excitations. The model has also been proposed to describe the edge states in 
fractional quantum Hall devices |lj . Tunneling through a point contact is then a practically 
ideal situation for comparing theory with experiments 0, 0, M. 

While most of the attention has focussed on transport properties so far, static prop- 
erties are also of interest. In particular, the effect of the impurity on the charge density, 
the equivalent of Friedel oscillations ||, has recently been considered 0. A very similar 
problem, in antiferromagnetic spin chains, was also discussed in J7J. See also for other 
related work. 

The ID Luttinger model with an impurity is integrable, and ultimately, using ad- 
vanced techniques of quantum field theory, ||, these Friedel oscillations should be exactly 
computable for all couplings, but there are important technical difficulties. The particular 
point g = 1/2, the equivalent of the "Toulouse point" in the anisotropic Kondo prob- 



lem, is equivalent to free fermions |T0[. It should thus be amenable by more elementary 
techniques. However, even in that case, exact expressions have not yet been obtained 
because the density operator is not local in the Fermion basis. We show in this paper how 
to circumvent this difficulty, largely based on a work of Chatterjee and Zamolodchikov 



2. Generalities 

Let us start from a description of the quantities involved in the Luttinger liquid, 



following closely @. Using standard bosonization formula [12 1 , one can write the electron 
creation operator in the spinless case - to which we restrict here - as a combination of two 
bosonic fields, (f>(x), 9(x). Decomposing into left and right moving parts, one has : 

ipix)^ resp. ip(x)^) oc exp [±ikfX + i\/7r(d(x) ± 0(x))] . (2-1) 

The canonical momentum for <p is d x 6 and they obey the commutation relations : 

[9(x),4>(x')] = ~ l -sgn(x-x'). (2.2) 
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The Friedel oscillations is describing the charge or density oscillations of the electrons in 
the presence of a barrier. The density operator [Tcj] : 

1 kf 

p(x) = p + ^d x (p + cos[2k f x + 2y/n^(x)] (2.3) 



with po = — the background charge. The hamiltonian for these fields is : 

/oo 
dx [U 2 + (d x (/>) 2 ]+H int . (2.4) 
-oo 

The term Hi n t describes a screened Coulomb interaction. The impurity will be coupled to 
the fields at one point x — by the term : 

H imp = ^d x (f)(0) + Acos[2V^0(O)]. (2.5) 

We will restrict to the case where Hi nt is short range, leading to a Luttinger liquid. In 
that case, the effect of the interaction is to renormalize the fields. The hamiltonian can 
then be brought into the usual form (setting vf — g) ■ 



fOO 

H = I dx 



oo 



87rsn 2 + -L(a^) 2 

87Tfir 



+ Acos[0(O)], (2.6) 



while the Friedel oscillations, subtracting the background charge density, are encoded into : 

< p{x) - po > = ^ kpx + < > ^ 

Po 

Here, the additional phase shift r]p = i n the cosine term arises from the unitary 

transformation that cancels the d x (j)(0) term in ( |2.5|) to get ( p.6|) . 

To proceed, we perform some manipulations. Decompose the field into its left and 
right components, <f> = 4>l + 4>R- Introduce then the left movers : 

(j) e (x + t) = -1= [<t> L (x, t) + <f> R {-x, £)] 

4>°(x + t) = —[ < f >L {x,t)-<t> R { r x,t)\ 

Now the impurity interaction reads Hi mp = A cos[v^2^ e (a; = 0)], while the observable we 
are studying is proportional to : 

cos <p(x, t) = cos ( —= [4> e + (x + t) + —j= [4> e - <f>°] (-X + t) 
lv2 V2 



Set now : 
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V2cj) e (x + t), x < § e R = V2(j) e (-x + t), x < 



$° = v/20°(x + t), x < = -V20°(-a; + 1), x < 

We then fold the system and recombine these left and right components into a single field 
to get the hamiltonian : 

H = H° + H e , (2.8) 



with 



and : 



H° = - 
2 



1 



dx 



dx 



8^(n°) 2 + J-(a x $°) 2 

Sng 



8ng(U e ) 2 + -^(d x $ e ) 2 
8ng 



$°(0) = 0, 
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+ Acos-$ e (0). 
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(2.9) 



(2.10) 



We will often refer to (|2.10|) as the boundary sine-Gordon hamiltonian. Using parity 
together with the decoupling of the odd and even fields we have : 



< p(x) — po > l^o/ \ l^p/ \ 

— = cos(2/c F x + rjF) < cos °(x) >< cos -$ (x) > . 

Po 2 2 



(2.11) 



We shall obtain results in the continuum limit. Then, since there is no boundary 
coupling for the odd field, one has MM : 



< cos °{x) >oc ( - 



(2.12) 



where a is a lattice coupling , x » a (a is defined by the propagator < $°(x)$°(0) >= 
— 4<?ln(x/a) for the bulk theory. The exact proportionality factor in ( |2.12| ) will be worked 
out later). Regarding the even field, things are more complicated. On general grounds, 
one expects : 

1 / a \ fl/2 r /t\ 1-9 

(2.13) 



\ fay/ 2 „ 

< COS (x) >OC - F 

2 \x/ 



a 



where the function F goes to one for large values of the argument, and vanishes linearly at 
small values of the argument (this latter property follows from the perturbative analysis 
in the regime g > 1/2). Continuum limit results will hold for x >> a and thus A << 1/a, 
while the product Aa (- 9 is kept constant. For simplicity we set a = 1 in what follows. 



3. Exact Priedel oscillations at g = 1/2 and T = 0. 

The point g = 1/2 corresponds to free fermions in the bulk, the analog of the Toulouse 
limit in the Kondo problem. To compute the universal function F exactly in that case, we 
will proceed in two steps: (i) We will show that cos ^- can be expressed as the product 
of spin operators in two decoupled massless Ising models, one of them having a boundary 
magnetic field, the other having fixed boundary conditions, (ii) These spin correlators in 
the boundary Ising models will be computed based on a method due to Chatterjee and 
Zamolodchikov. 

In the following we do computations for the even field (ie the field having boundary 
interaction); results for the odd field follow simply by taking the limit of Dirichlet boundary 
conditions in the former solution. 

A natural way to relate the boundary sine-Gordon model to two copies of the Ising 
model is by using bosonization. At the free fermion point this was considered in [|15| , where 
it was shown that the boundary interaction Acos^$ e (0) becomes linear in the complex 
Dirac fermion fields, ijj± : 'ip±. (In the conformal limit, ijj± and ^ ± are the left and right 
components and ± is the U(l) charge.) The boundary equations of motion at x = take 



the form 15 



■?/>+ + if>- - i>_ - ^+ = o 

id t {tf>- - ^ + ) - A 2 $L - <M = (3.1) 
zd t (V>+-?_) + A 2 (^ + -^ + ) = 0, 

where A oc A and A is defined in ( p.lOQ . Note that these equations are valid regardless of 
whether there is a bulk mass term or not. 

Defining real components of the fermi fields as ijj± = ± iif)2, ij)± = —(ipi =t ^2)5 



one finds that (|3.1] ) is equivalent to 



(3-2) 

za t (^2-^ 2 ) + A 2 (^2 + ^ 2 ) = 0. 



One can now compare these equations with the boundary Ising equations of motion |TT[ 



one sees that at the free fermion point, the boundary sine-Gordon theory is equivalent to 
two boundary Ising models, one with fixed boundary conditions (equivalently a magnetic 
field h = 00), the other with a varying magnetic field h oc A. The precise relation between 
h and A is however more delicate to obtain, see below. 
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A less natural, but more precise way to relate boundary sine-Gordon model to two 
copies of the Ising model is to use a scattering description, based on the boundary S-matrix 
and form factors. To clarify some of the following discussion, let us for a while add to H e a 
term proportional to f_ cos$ e (x), generalizing it into a massive boundary sine-Gordon 



model, which is well known to be integrable [^6[. At g = 1/2, to which we restrict, the 



spectra of the theory is composed of solitons and anti-solitons (there is no breathers) and 
the boundary reflexion matrix is described by the two terms, P and Q given by [TO : 



R± = P, P±= Q, (3-3) 

where the labels ± refer to the soliton and antisoliton of the sine-Gordon spectrum. For 
our purposes it is much more convenient to take the following combinations : 

P + Q = i tanh ( i 

. / 7T u\K — i smh t) 
P — Q = ztanh i- 



4 2 J k + i sinh 9 

The key remark here is that the combination given are just reflection matrices for the 
fermion of the Ising model in a low temperature phase in a boundary magnetic field h, 



with k = 1— 4?Th . (ft = r in notations of [0). One combination describes fixed Ising spin 



boundary conditions (or infinite boundary magnetic field), and the other is at finite Ising 
boundary magnetic field. This suggests that in the basis of symmetric and antisymmetric 
combinations of solitons and anti-solitons, the theory decouples to two Ising models with 
boundary magnetic fields. 

This can be confirmed for instance by considering form factors, ie matrix elements of 
the operators in the quasiparticle basis |p.7|| . The form factors for the operators e ±J */ 2 in 
the sine-Gordon model at g = 1/2 are given by : 

2n 

f ± (e 1 ,...e 2n ) ei _ e2n = (±2i) n ^ir^l[uf^ 2 l[(e l u l -e J u j y^ (3.5) 

i=l i<j 

with the notation Ui = e 6 \ 6 the usual quasi particle rapidity |T7|j . In order to show the 
relation to the Ising model, we change the basis from \A>, \S > to -^(\S > +\A >) and 
-^(|>S > —\A >). The explicit calculation is done in the appendix, in the new basis the 
form factors are : 

4..., eL (0i, 2n ) = i n (-ir~/ 2 Jltanh (^=^) ' ^ J/ , (3-6) 

i<j ^ ' 
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where the isotopic indices refer now to the sign of the combination between solitons and 
anti-solitons and n_ denotes the number of indices having value — . This formula holds 
only for n_ even, which will be the case for the correlator of interest. Recall now the 



expression of form factors for spin operator a in the Ising model |18| 



/"((?!,..., e 2n ) = i n H tanh ~*~2~^> (3 - 7) 

the form factors of (jl on an even number of particles vanishing. 

In the presence of the boundary, the one point function of of interest is obtained 
as < B\ cos - i 3 ^ |0 >, where \B > is the boundary state, |0 > the ground state. The 



boundary state can be expressed in terms of multiparticle states following [16] . In the new 



basis mixing solitons and anti-solitons, this boundary state actually factorizes : 

\B> = \B >+ ®|J3 >_, 
\B >± = exp f / ^A* ± (-(3)A* ± (f3)K ± ((3) 



(3.8) 



where A*, denote the creation operators in the new basis, and K±{6) = ±(P±Q)(^ — 9). 
Here, the additional ± sign occurs because the K matrix is obtained through the general 
formula |L6| K ab (9) =R b -( l f - 9); soliton and antisoliton are conjugate in the sine-Gordon 



model, while the Ising fermion is self conjugate. Expanding these boundary states, using 
the form factors ( |3.6| ) and ( |3.7| ), we see that the correlation function becomes indeed a 
product of two spin operators (in the low temperature phase of the Ising models, where 
the one point function of the disorder operator vanishes) : 

< cos — >=< a > h < a >oo . (3.9) 

Having shown (|3.9| ), we may now let the mass of the sine-Gordon model (the amplitude 
of the bulk cos $ e term) go to zero to recover the original problem. The correspondence 
between the variable A in the boundary sine-Gordon, and the field h in the boundary Ising 



model, can be found in that case using results of |n|. Using the conventions of |TT| for 
the boundary Ising model one finds Tg = 8nh 2 = 4ir\ 2 so h = This completes our 
derivation of the relation between the Friedel oscillations and the Ising model. We can 
now concentrate on finding the correlation of the spin operator in the Ising model with a 
boundary magnetic field. 
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At T = 0, the one point function of the spin with a boundary has been evaluated 
in the very interesting paper JO]. There it is shown that (we trade the variable A of the 
original action for h), introducing the variable X = 4nh 2 x, and setting a(X) = {a(x)) h , 
the following holds : 



d 2 



dX 2 



+ 1 Jc 



d / 1 9 1 
dX + \X + 16X2 



a(X) = 0. 



The solution of this equation that describes the appropriate physics is : 



(a(x)) h = 2 13 / 8 v^F/ix 3/8 ^(l/2,l;87r^ 



x 



(3.10) 



(3.11) 



where a is normalized as usual in the bulk. Here, $ is a degenerate hypergeometric 
function. It is simply expressed in terms of Bessel functions as \I/(l/2, 1; 2x) = -^Ko(x). 
At short distance one finds that a rs — 2 13 / 8 /i x 3 ^ 8 \nx and at large distances o ~ 2 1 / 8 x -1 / 8 . 
Taking into account the fact that the physical observable in the case of Friedel oscillations 
is ( [2. 11|) one finds : 



< p(x 



5 — = A^h cos(2k F x + r7 F )*(l/2, 1; 8nh 2 x), h = 4=- 
Po V2 



The asymptotic behaviour is ^*>-<«>' oc — lnx cos(2/cfx+?7f) as x 



<p(x)-p > 
PC 

x~ x l 2 cos(2kFX + t]f) as x — > oo. 

Our result agrees with the numerical simulations in 



(3.12) 
Oand <gOgWo> 



oc 



except at very small x, where 



it was found that ^ P{X> po Po ^ behaves as a power law, with a small negative exponent. This 



<p(x)-p > 
PO 

discrepancy is very likely due to fact that, for stability reasons, the true UV region is 
difficult to access numerically - indeed, at intermediate values of x, our result does behave 



like a power law [20 



4. Exact Friedel oscillations at g = 1/2 and T ^ 0. 



We can also extend the computation of |Oj to finite temperature. The main idea is 
still that the boundary magnetic field does not destroy the free field structure of the Ising 
model ||21| , leading to a determination of the correlator by elementary considerations |Tl 



For convenience, we first rotate the geometry so now the boundary lies along the x 
axis. As mentioned earlier ( |3.1| ), the equations of motion for the Majorana fermion of the 
Ising model with boundary magnetic field are (in imaginary time) : 

[{d z + 4tnh 2 )ifj(z) - {d- z - 4inh 2 )ij(z)} z= - z = 0, (4.1) 
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where z = x + iy. Having this condition, the idea is to introduce a fermion field x( z ) = 
(d z + 4:iirh 2 )ifj(z) and its conjugate and realize that the previous boundary condition 
states that x(z) is the analytic continuation of x{ z ) to the lower half plane. Hence the 
correlator < x(z)/i(u>, u>) > is analytic in the full z plane with two square root branch 
points at z = w and z = w. 

With a finite temperature, the argument is the same apart from the periodicity in 
the imaginary time direction, which after having rotated the system is the x direction, is 
z — > z + -g. The boundary condition, a local statement, remains the same. 

The next step is to write a global form for the correlator < x( z )l^( w 7 w) >, and this 
is where the effect of temperature will be seen. On the cylinder, one requires the right 
hand side to have square root branch points at w + nT, w + mT, n, m integers, and to be 
periodic. One can therefore write : 

... s /sin7rT(z — w) \ /sin tvT(z — w) x ' "~ 



\ ttT J \ tvT 

= tvT cot[TvT(z — w)] A(w, id, T) + tvT cot[-7rT(z — w)] A(w, id, T) + B(w, id, T). 

(4.2) 

Observe how the right hand side is periodic in z — > z + T. In this form, the coefficients, 
A, A, B are unknown that need to be fixed. We can now use the operator product expan- 
sions, for example in (z — w) (uj = e l7r / 4 ) : 

x(z)fJ,(w, w) = —7={z — w)~ 3 / 2 { — (j(w, w) + (z — w)(2d w + Amh 2 )a(w, w) 

V2 2 (43) 

+(z - w) 2 (4dl + 16iTvh 2 d w )a(w, w) + . . .}. 



and a similar relation for x- The reader might fear that ( [4.3|) which holds in the plane 
might be changed by (3 dependent terms when the geometry is compactified. Actually, it is 
well known PI2| that short distance expansions are invariant in conformal mappings. This 
constraint for instance determines the coefficients 4, 8/3 in the operator product JTJJ : 



ip(z)fj,(w,w) = -j=(z-w) 1/2 



a(w, it?) + 4(z — w)d w a(w, w) + —(z — w) 2 d1 ] a{w, w) + . . . 

(4.4) 



Now taking the global form (|4.2|) and expanding in (z — w) we get relations between the 
coefficients A, A, B and the spin function and its derivatives (similarly expanding in (z — w) 
and comparing with the OPE of x/i we get more relations). Then, by matching ( f4.2| ) and 



and eliminating the unknown coefficients one finds that the one point function of the 
spin satisfies the equation : 

— (7rT) 2 (coth27rTy) 2 coth27rTy 



+ 



< cr(y) >= 0. 



In particular, this equation reproduces eq. (30) of [[TI| in the limit T —>■ 0. To solve this 
equation, let us set < a{y) >= [sinh(27rTy)] ' f(y). Introduce then the new variables 
A = ^h- a nd Y = (1 — coth27rTy)/2, and set f(y) = f(Y). We then find an hypergeometric 
equation : 

( y - r2 >(W + ( 1 + I- 2y )5F-l} /(y) ^ ' (4 ' 6) 

The physical solution of the problem , going back to the original x variable, is then : 

( AnT \ 1/8 /l 1 h 2 l-coth27rTx\ 
< ^ )> ^(sinh2^J F U'2 ;1 + V 2 J' (47) 

To select this solution, we observe that the one point function at large x can should 
expand, by conformal invariance, as a sum of exponential terms exp(— 4nTxA), with A 
conformal weights of the (central charge) c = \ conformal field theory. This excludes the 
other hypergeometric function solution of (|4.6|) . The normalisations are fixed by looking 
at the case h — > oo. If h — > oo, one has fixed boundary conditions for the spin. The 
one point function can then be obtained by conformal mapping from the half plane where 
it is known to be []S3[ < a(x) > hp = 2 1 / 8 x -1 / 8 (with < a(x)a(0) >= x -1 / 4 in the full 
plane). Setting z' = l n fxj maps the half plane on the semi-infinite cylinder, the line 
boundary Rez = becoming the circle boundary Rez' = 0,lmz' G [— 1/2T, 1/2T]. In 

1/8 

this transformation one finds for the semi-infinite cylinder < cr(x) > s ic= ( ^~Tnr~T~ ) > ^ n 
agreement with ( [4.7| ). Also, as x —>■ oo, we expect to recover the same result since at large 
distance the spin sees fixed boundary conditions. This fixes the normalization constant in 



Finally, we can study the limit T — > 0. By using the standard transformation formulas 
z — > l/z for the argument of the hypergeometric functions, together with the definition of 
\I/ in terms of the basic degenerate hypergeometric function $ (all notations are those of 
PHI), one finds, as 7, z — » 00, j/z finite 

F(a,a;7,*)« (^) #(1/2, 1; - 7 /z), 



where we used 

F(a, a + 1 - 7; 1, 1/z) ss $(a, 1; -7 A), 

which can be proved using the series representation of the involved functions. Therefore, 
as T — > 0, one finds < a(x) 2 1 / 8 x~ 1 / 8 (87r/i 2 x) 1/ ' 2 \l/(l/2, 1; 8nh 2 x), in agreement with 
the foregoing results [JTI]. 

The qualitative effect of the temperature can be seen on the small and large x limits 
of the one point function. One has 

<a(x) >^-2 9/8 r 1/a r ^ + ^ , x«l/T (4g) 
< > h « (87r) 1/8 e-" Tx/4 , x » 1/T. 

We can now come back to our original problem, the Friedel oscillations. Using ( |2.11j ) 
with the previous solution we find the final expression : 



<p(x)-po> , n , ,f 4nT \ 1/2 /l 1 „ h 2 l-coth27rTx\ 

(4.9) 

where as before h = The periodicity x — > x — i/T of the final solution might appear 
a bit surprising since it is the imaginary time y which is compactified on a circle of radius 
1/T, not x. It can be understood from the fact that one is dealing with a massless theory 
in the bulk. In fact, this periodicity can also be seen by taking the massless limit of the 
expressions for finite temperature correlators obtained in ||25|| by using the form factors. 

We finally observe that the same technique could be applied to study the screening 
cloud at the Toulouse point of the anisotropic Kondo problem |26| . This will be described 
elsewhere. 

Acknowledgements: We thank R. Egger for pointing out this problem to us, for 
many useful conversations, and for kindly checking our formula against his numerical 
results. 



Appendix A. Change of basis for the sine-Gordon form factors. 

In this appendix we want to show how the one particle form factor for the operator 
e ±*<^/2 j n sine-Gordon model at the free fermion point is related to spin form-factor in 
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the Ising model. The form factor is a tensor function in the space of isotopic indices and in 
the new basis described by -^(\S > +\A >) and -^(\S > —\A >) (which we will denote 
by + and — in the following), the form factor is simply given by : 

/ ell ... l4n (01, M = All ' ' ' A £Al,-^n(01, -> *2»), (A.l) 

where the primed indices denote the new basis and with : 

_ 1 /l 1 



a '-=tAi -i) (A - 2) 

describing the change of basis. The forms factor in the \S >, \A > basis was given before : 

2n 

/ ± (e l ,-^2n)6 1) ... ) e an = {±2l) n ^/ Ul ...U 2n J] U^ 2 - '/"/)' , ' / (A-3) 

where the non-zero form factors are those preserving charge and we used the variables 

Ui = e 9i . 

Let us first look at the form factor with all primed indices chosen to be +, it is given 

by : 



/i.. + (01,-.,02n)=(±t)' 



2n / \ A\ 

By making a dilatation one sees hat this is an homogeneous function of degree zero. This 
function has poles for all values of itj = — Uj and zeroes for Ui = uj. Only the sign has to 
be worked out, the final form factor is : 

/£..+ (*!, 0*0 =i n Ht^h( d -^A . (A.5) 

Imagine now that some of the primed indices are — , then the function is still homegeneous 
of degree zero but if let say <pi are the rapidities of the particles having indices — , it is not 
difficult to show that there are now no zeroes and poles between terms involving 9i and 
<pj. There are still zeroes and poles between each type of particle and we find that the 
form factor is given by : 

/a — a \ K+ e j-l/ 2 

/± ^(0!,...,^) =<»(-l)»-/ a J]tanhf ^Y^) ■ (A - 6) 

i<j ^ ' 

In this expresion, ri- is the number of indices e[ of type — ( ([A.6 ) holds only when n_ is 
even, which is the case for our computation). 
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